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Abstract 

A study is made of charm and bottom flavor matching conditions on 
parton densities. Starting from an MS three- flavor density set, where the 
scale /i < rUc, and using the recently derived two-loop matching conditions, 
we provide a new set of four-flavor parton densities where rUc < n < ruh and 
five-flavor densities where fi > mf,. The effect of the next-to-next-to-leading 
order matching conditions on the evolution equations is important for scales 
just above the transition regions. This includes the small x and small 
domain studied by the HI and ZEUS experiments at HERA. At small x the 
effects of the matching conditions never die away even for large fi"^. 

PACS numbers: ll.lOJj, 12.38Bx, 13.60Hb, 13.87Ce. 
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1 Introduction 



Quantum chromo dynamic predictions for experimental cross sections and 
distributions in perturbation theory rely heavily on accurate knowledge of 
parton densities. Several groups [|l|, 0, have extracted these densities 
from global fits to data with the latest theoretical information on MS coeffi- 
cient functions. At present the evolution of these densities via the Altarelli- 
Parisi (AP) equations 0] uses the information on leading order (LO) and 
next-to-leading order (NLO) splitting functions. Unfortunately the next- 
to-next-to-leading order (NNLO) splitting functions are not known but some 
interesting pieces of information are available ||^, The complete split- 
ting functions should be known soon. 

The description of heavy quarks within this analysis recently received a 
lot of attention due to the data on deep inelastic production of D* mesons 
from HERA [l^. The global fitting groups have adopted different ap- 



proaches. The CTEQ5 analysis describes charm and bottom densities via 
the so-called ACOT prescription |Tl|, which is a one-loop matching condi- 



tion between three-flavor and four-flavor densities at the scale /z = ttIc and a 
corresponding one-loop matching between the four-flavor and flve-flavor den- 
sities at the scale /i = m^. This is not done in the MRST density sets. Instead 
they impose matching conditions that the logarithmic derivative of the deep 
inelastic structure functions with regard to the scale should be continuous, 
see |12| for details. This yields different charm and bottom densities. The 



GRV group flSl adopt the approach that one does not need any densities 
other than a three-flavor MS set because the convolution of these densities 
with the NLO heavy quark coefficient functions provided in yields an 
excellent fit to the presently available data on F2^c{x-, Q^-, 'mc)^ which is stable 
under scale variations. By never taking the limit that rric —>■ the theoret- 
ical prediction has no collinear singularity problem. Note that the physical 
threshold for a heavy quark antiquark pair is at Q'^{1 —x)/x = Am?. However 
the physical threshold is distinct from the matching scale where one switches 
between parton density sets. 

The parton densities with nj and + 1 fiavors are related by a set 
of operator matrix elements (OME's). The order a1 OME's were recently 
derived in They contain terms with ln*(/i^/m^) i = 1, 2 as well as non- 
logarithmic terms. They have the property that the n/ + 1 fiavor densities 
vanish in LO and in NLO when the scale n = m. In NNLO there are finite x- 
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dependent discontinuities at this scale, which we refer to as NNLO matching 
conditions. In this respect the matching conditions on the parton densities 
are similar to those for the two-loop running coupling constant derived in 



1^, |2y]. These NNLO matching conditions can be important numerically. 
Note that the running coupling constant is not small at the scale /i = m^. 
The charm density constructed from the two-loop matching conditions has 
been used in a recent study of variable flavor number schemes for the charm 
component of the deep inelastic structure functions [1^. Here we would hke 
to present a complete set of parton densities for light (u,d,s,g) and heavy 
(c,b) partons which satisfy the NNLO boundary conditions and are evolved 
with NLO splitting functions. When the three-loop anomalous dimensions 
become available they can be included in this analysis. At present our study 
primarily focuses on the effects of the discontinuous matching conditions 
and the differences between the resulting three, four and five-flavor parton 
densities as the scale /i increases through the regions < /i < rrih and 
fi > mf, respectively. 

In Sec. II we present some technical details about our choice of input 
parton densities and the two-loop matching conditions. All densities are 
derived in the MS scheme. Then in Sec. Ill we give plots of the parton 
densities. Finally comparisons are made between the presently available NLO 
sets and our (NNLO) set. A discussion of the solutions of the evolution 
equations is given in the Appendix. We intend to make our computer package 
for these densities available in due course. 
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2 The parton densities 



We present a consistent set of MS parton densities containing three, four 
and five flavors for scales satisfying fi < rric, rric < rnf, and fi > rrib 

respectively. The evolution of the densities is done with our own computer 
code written in C++ and some details are given in the Appendix. The code 
uses the direct x-space method to solve the evolution equation M similar to 



that in and allows us to evolve both light and heavy parton densities 
in LO, NLO and NNLO (the latter using the NLO weights). The weights 
for the calculation are computed analytically from the LO and NLO [§] MS 
splitting functions thus removing possible instabilities in the numerical inte- 
grations. Hence the program is very efficient and fast. The results from the 
evolution code have been thoroughly checked against the tables in the HERA 
report |18[. We use weights in LO and NLO for rif = 3, 4, 5 for evolving the 
gluon fg{nf,x,fi'^), the singlet quark densities fg^UfjXjfi"^), the non-singlet 
valence quark densities fk-k{nf, x, iJ?) and the non-singlet sea quark densities 
f^^{nf, X, /i^). As the scale increases across the charm and bottom matching 
points the sets are redefined to include densities for the c and b heavy quarks. 
The program allows us to use LO, NLO or NNLO matching conditions for 
the generation of the heavy fiavor densities. 
The number densities are defined as 

fk-kinj, X, /i^) = fk{nf, X, /i^) - /fc(n/, x, ^?),k = u,d 
fk+ki^f, X, /i^) = /fc(n/, X, /i^) + /fc(n/, x, iJ^),k = u, d, s, c, b 

(2.1) 

f^{nf, X, /i^) = fk+kijif, X, /i^) (2.2) 

k=l 

fg^irif, X, ^^) = fk+-k{nf, X, /i^) - —f^{nf, x, /i^) . (2.3) 

We start our LO evolution using the following input from at /Xq = 
ftlo = 0.26 (GeV/c^)2 

u(3, X, /ig) = X1t^(x, /J-lq) 

= 1.239 x°-^^ (1 - xf-^^ (1 - 1.8y/x + 9.5x) 
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= 0.614 (1 - xu^{x, filo) 

x{fd{^, X, nl) - /ii(3, X, /io)) = xA{x, /iLo) 

= 0.23 (1 - xY^-^ (1 - 12.0v^ + 50.9x) 

a;(/j(3,x,/io) + /s(3,x,/io)) = x{u + d){x,filo) 

= 1.52 x°-^5(l-a;)^-^ (l-3.6v^ + 7.8a;) 

a^/g(3, a^, /^o) ~ xg{x,fii^Q) 

= 17A7 x^-^ {1-xf-^ 

a^/s(3, a;, /Xq) = x/g(3, X, /ig) = xs(x,/iLQ) 

= xs(x,/i2o) = 0. (2.4) 

Here A = d — m is used to construct the non-singlet combination. We start 
the corresponding NLO evolution using the following input from at /ig = 
/x^LO = 0.40 (GeV/c')2 

X fu—u(3: X , ^q) = XUy(^X, fij^j^o) 

= 0.632 (1 - a;)^-°^ (1 + 18.2x) 
a^/(i-(j(3, a;, /ig) = /i^Lo) 

= 0.624 (1 — xY''^ xUy{x, AiNLo) 

a;(/j(3, X, /ig) - /ii(3, x, /ig)) = xA(a;, /i^Lo) 

= 0.20 (1 - a;)i2.4 _ 13.3^ + gO.Ox) 

a;(/<j(3,x,/io) + /s(3,x,/io)) = + J)(x,/iNLo) 

= 1.24 x°-20(l -x)*-^(l -2.3v^ + 5.7a;) 

a^/g(3, a;, /ig) = a;5f(x, /x^lo) 

= 20.80 x^-^ (1 - a;)^-^ 

a^/s(3, X, /ig) = x/g(3, X, /ig) = xs(x, /i|v^LQ) 

= xs(x,/i^Lo) = 0. (2.5) 

From the above densities we form the combinations that we evolve and step 
across thresholds. The NNLO densities for /i^LO < 1^^ < f^^l are replaced with 
NLO densities. The heavy quark masses rric = 1.4 GeV/c^, mi, = 4.5 GeV/c^ 
are used throughout the calculation together with the exact expression for the 
running coupling constant a;s(/i^), represented as the solution of the following 
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differential equation 



or in the imphcit form 



In 



47r 



/3i 



(2.6) 



(A 



("/) ^|2 

EXACT 



_ ^ 

)2 /?oas(/i2) /3o 



In 



47r ^ A 



(2.7) 



where /3o = ll — 2r2//3 and /9i = 102 — 38r;,//3. The values for A^^^ are 
carefully chosen to obtain accurate matching at the scales ml and respec- 
tively. We used the values hitler = 299.4, 246, 167.7, 67.8 MeV/c^ respec- 



tively in the exact formula (which yields a 



EXACT/ 



mr 



0.114, «!^XACT, 



mt 



0.205, a 



EXACT/ 



ml) = 0.319, a^ACT^^ 



NLO/I 



0.578 



and Ag^'^-') = 204, 



0.125, 



mt 



0.232, aLO(m2) = 0.362, aLO(/i2^) = 0.763 ) for the LO formula. Note that 



175, 132, 66.5 MeV/c respectively (which yields a\ 

[ml) = 0.362, a^^if^lo) 
we have not used the two-loop matching of the running coupling constant 
asiiJ''^) at the same scales from [^,|^ to focus on the matching conditions 
on the flavor densities. Numerically the discontinuity in the running cou- 
pling constant across the charm threshold is aproximately two parts in one 
thousand, which is far too small to affect our results. 



Three flavor evolution proceeds from the initial /Xq to the scale /i^ 



mz 



1.96 (GeV/c ) . At this point the charm density is then defined by 

U,{nf + 1, ml) = aliuf, m^) lAHil) ® /^{nj, ml) 



+A 



Qg\ 



fginf,ml) 



(2.8) 



with Uf = 3 and = as/ An. We have suppressed the x dependence to 
make the notation more compact. The ® symbol denotes the convolution 
integral /^f? = / f {x/y)g {y)dy/y, where x < y < 1 The OME's A^^^fi"^/ ml), 
Aqglfi"^ /ml) are given in [|1^. The reason for choosing the matching scale 
fi at the mass of the charm quark is that all the \n{fi'^/ml) terms in the 
OME's vanish at this point leaving only the nonlogarithmic pieces in the 
order al OME's to contribute to the right-hand-side of Eq.(2.8). Hence the 
LO and NLO charm densities vanish at the scale fi = mc- The NNLO charm 
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density starts off with a finite x-dependent shape in order a^. Note that we 
then order the terms on the right-hand-side of Eq.(2.8) so that the result 
contains a product of NLO OME's and LO parton densities. The resuh is 
then of order and should be multiplied by order coefficient functions 
when forming the deep inelastic structure functions. 

The four-flavor gluon density is also generated at the matching point in 
the same way. At /i = rric we define 



fUrif + 1, 



+al{nf,ml 



+A 



gq,Q\ 

fgi^f^ 



l)®f(nf,ml 



(2.9) 



The OME's 



(/i^/m^), A^gQ^fi'^/ml) are given in |T5|. The four-flavor 
light quark (u,d,s) densities are generated using 



+aJn 



iNS 



l)®h+dnf,ml). (2.10) 



The OME ^^^^^(/iV^c) is given in fTSlI and the total four-flavor singlet quark 



density in Eq.(2.2) follows from the sum of Eqs.(2.8) and (2.10). The non- 
singlet density then follows from Eq.(2.1). In Eqs.(2.9) and (2.10) nf = 3. 
The remarks after Eq.(2.8) are relevant here too. 

Next the resulting four-flavor densities are evolved using the four-flavor 



weights in either LO or NLO up to the scale n = 20.25 (GeV/c ) . The 
bottom quark density is then generated at this point using 



fh+-b{nf + l,ml) 



l(S) 



(2.11) 



and the gluon and light quark densities (which now include charm) are gen- 
erated using Eqs. (|2.8| )-( |2.10| ) with nf = A and replacing by m^. There- 
fore only the nonlogarithmic terms in the order OME's contribute to the 
matching conditions on the bottom quark density. Then all the densities are 
evolved up to higher fi"^ as a five-flavor set with either LO or NLO splitting 
functions. This is valid until fi = nit ^ 175 GeV/c^ above which one should 
switch to a six- flavor set. We do not implement this step because the top 
quark density would be extremely small. 
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The procedure outlined above generates a full set of parton densities 
(gluon, singlet, non-singlet light and heavy quark densities,) for any x and 
/x^ from the three-flavor LO and NLO inputs in Eqs.( ^.4D and ( |2.5| ). Note 
that one could also use the formulae above in fixed order perturbation theory. 
In this case the four-flavour densities are defined by extending the integrals 
on the right-hand sides of Eqs.(2.8)-(2.11) to 




(2.12) 



and 



+a,(n/,/i 



mt 



fk+ki'^f + '^'f^^) = fk+ki^f^f^^) 



2\ 4NS / 



fk+kirif,ix^) 



(2.13) 



(2.14) 



for rif = 3 and ^ /^^ < Then the five-flavor densities are deflned by 



Till. 



I flu 
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(2.15) 



for Uf = 4 and > m^. Also one should replace n/ = 4 and by ml in 
Eqs.(2.12)-(2.14). In this case no four-flavor or flve-flavour evolution is made 
so the logarithmic terms in /x^/m^ and/or n^/ml are not summed. We will 
show the differences between this fixed order perturbation theory (FOPT) 
treatment and the evolved treatment in the next Section. 
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3 Results 



Here we present results from the evolution of the parton densities. The inputs 
are the three-flavor densities at /iq in Eqs.(2.4),(2.5) which are evolved up to 
the scale = rric = 1.4 GeV/c^. During this evolution the number of light 
flavors n/ = 3 in both the MS splitting functions and the running coupling 
constant. 

We start by giving the four-flavor densities, where = 4, in the region 
between rric < < frif, which follow by evolution from the matching con- 
ditions in Eqs. (2.8)-(2.10). We present results at the scales /x^ = 1.96, 2, 
3, 4, 5, 10 and 20 in units of (GeV/c^)^. First we show the charm density 
a;c^^'"°(4, X, /i^) in two ranges (a) 10"^ < a; < 1 and (b) 10~^ < a; < 1 in 
Figs. 1(a) and 1(b) respectively. We notice that this density starts off nega- 
tive at small x but it is positive at large x so that the momentum sum rule 
is satisfied. To show the effect of resumming the logarithmic terms via the 
evolution equation from the charm threshold as compared with just comput- 
ing the integrals in Eq. (2.12) at all scales fi we show in Fig. 1(c) the ratios 
R^^^^{x, fi"^) = xcevolved(4, x, /i^)/xcFOPT(4, X, /i^). Here FOPT stands for 
fixed order perturbation theory. The effects of the evolution are especially 
significant at small x and large x. Notice that the discontinuity at x ~ 0.01 
is caused by the change in sign of the NNLO charm quark density. For a 
comparison we have also shown the NLO results from the MRST98 set 1 
and CTEQ5HQ [|l|] parton density sets in Figs. 1(d) and (e). These groups 
use different input densities so a direct comparison does not have any true 
significance. Nevertheless our density is larger than the MRST98 set 1 result 
and smaller than the CTEQ5HQ result at small x and large n^. 

In Fig. 2(a) we show the four-flavor gluon density x(y'^'^^'^(4, x, /i^) in the 
same range 10^^ < x < 1 for the same scales as in Fig.l. We also show in 
Fig.2(b) the ratios i?J^^°(x, /i^) = x5(evolved(4, x, /x^)/x5fFOPT(4, z^^) for 
the same scales, where we use Eq.(2.13) for the FOPT density. The effect of 
the suppression of the charm density at small x translates into an increase 
of the gluon density at small x. For comparison we show the three-flavor 
NLO gluon density in Fig. 2(c). We have also shown the NLO results from 
the MRST98 set 1 and CTEQ5HQ parton density sets in Figs. 2(d) and 
(e). These densities do not increase as rapidly at large yU^ because they use 
different inputs. 

In Fig. 3(a) we show the singlet quark density xS^^'"°(4, x, /x^) in the 
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range 10^^ < a; < 1 for the same scales as above. Then is Fig. 3(b) we show 
the ratios R^^^'^{x, n'^) = a;SEvoLVED(4, a;, /x^) /a;SFOPT(4, a;, /x^), where we 
use Eq.(2.14) for the FOPT density. This ratio shows increases or decreases 
depending on the x and /i^ values. It is appreciably smaller at large /i^, 
which reflects the differences between the three-flavor and four-flavor gluon 
densities. For comparison we show the three- flavor NLO density in Fig.3(c). 

Next in Fig.4(a) we show the nonsinglet density xa^^^'^{4:,x, n'^), where 
a — {u + u)/2, in the range lO""^ < x < 1 for the same scales. In Fig. 4(b) 
we show the ratios E!^^^'^{x, ii^) = a;crEvoLVED(4, x, /i^)/x(TFOPT(4, a;, ^^) for 
the same scales, where we use Eq.(2.14) for the FOPT density. The ratio 
is significantly below unity at large x. In Fig. 4(c) we show the three-flavor 
NLO result. The difference is small. 

We complete our presentation of four- flavor densities by showing in Fig. 5 (a) 
the strange quark density a;s^^'"°(4, x, /i^) in the range 10~^ < x < 1. In 
Fig.5(b) we show the ratios i^^NLO^^^ _ ^^^evolved(4^ ^2) /xs^^^^{A, x, /z^) 

for the same scales. In Fig. 5(c) we show the three- flavor NLO density, where 
again the difference is small. We have checked that these densities satisfy 
the momentum sum rule for four flavors. 

Now we move up in scale to consider /i > rrib — 4.5 GeV/c^ which is the 
five-flavor region. The parton densities in this region are now generated from 
the previous four-flavor set by using the conditions in the Eqs.(2.8) -(2.10) 
with n/ = 4 at /^^ = ml and replacing ml by mf. Here we show plots for the 
scales iJ,^ =20.25, 25, 30, 40, and 100 in units of {GeV/c^f. 

The first density to consider is the bottom quark density. We show in 
Fig. 6(a) xb^^^'~'{5,x, fi'^) in the range 10~^ < x < 1 for the scales mentioned 
above. Notice that it is negative for small x and small /x^. The region 10"^ < 
a; < 1 is shown in Fig. 6(b) to demonstrate that the density is positive for 
large x. In this respect it is like the charm density in the region just above the 
four-flavor matching point. In Fig. 6(c) we show the ratios R^^^'^{x, /i^) — 
xb^^^^'^^^{5.,x,fi^)/xb^^^^{5,x,fi'^), where Eq.(2.15) is used for the FOPT 
density. Here it is clear that the effect of the evolution is appreciable at 
small X. For a comparison we have also shown the NLO results from the 
MRST98 set 1 and CTEQ5HQ parton density sets in Figs. 1(d) and (e). 
Note that these groups use different input densities so a direct comparison 
does not have any true significance. Nevertheless both bottom densities are 
larger than ours at small x and large /i^. 

The five-flavor charm density a;c'^^^'~*(5, a;, /x^) is shown in Fig. 7(a) in the 
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range 10^^ < a; < 1 for the same scales. In this case the ratios of the evolved 
densities to the FOPT densities is very close to unity for all x and /x^ values 
so we do not show a plot. For a comparison we have also shown the NLO 
results from the MRST98 set 1 and CTEQ5HQ parton density sets in Figs. 
1(b) and (c). The former has a smaller density at large //^ and the latter a 
larger density. 

Next we show in Fig. 8(a) the gluon density xg^^^'^{b,x, ii'^) in the range 
10~^ < X < 1 for the same scales. In this case the ratios of the evolved 
densities to the FOPT densities is very close to unity for all x and //^ values 
so we do not show a plot. In Fig. 8(b) we show the corresponding three-flavor 
NLO density. The latter is larger at small x and large /i^. 

In Fig. 9(a) we show the nonsinglet quark density a;S^^^*-'(5, a;, /x^), where 
a — {u + u)/2, in the range 10~^ < x < 1 for the same scales. Also in this 
case the ratios of the evolved densities to the FOPT densities is very close 
to unity for all x and values so we do not show any plot. In Fig.9(b) we 
show the corresponding three-flavor NLO density. 

In Fig. 10 (a) we show the nonsinglet quark density a;(j^^^°(5, x, ii'^) in the 
range 10"^ < x < 1 for the same scales. Also in this case the ratios of the 
evolved densities to the FOPT densities is very close to unity for all x and jj? 
values so we do not show any plot. In Fig. 10(b) we show the corresponding 
three-flavor NLO density. 

In Fig.ll(a) we show the strange quark density xs^^^'^ib^x, fi^) in the 
range 10~^ < x < 1 for the same scales. Also in this case the ratios of the 
evolved densities to the FOPT densities is very close to unity for all x and jj? 
values so we do not show any plot. In Fig. 11(b) we show the corresponding 
three-flavor NLO density. We have checked that the five-flavour densities 
satisfy the momentum sum rule. 

The above plots demonstrate that the NNLO matching conditions do in- 
fluence the parton densities appreciably in regions just above /i — rric and 
fi = mil. This has consequences for the analysis of HERA experiments be- 
cause a lot of the data is at small x and small values of the scale Q^. In fact 
all of the data for x < 10^"^ has Q"^ < 100 (GeV/c^)^. Even for the scale //. 
much larger than mf, the boundary conditions are still important. For // > 10 
GeV/c^ for example the rapid rise of the five-flavour gluon density means that 
it dominates over all the other parton densities at small x. There is roughly 
a ten percent difference between the three-flavor and flve-flavor densities at 
small X and large /x^, which can be important for precision phenomenology. 
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Appendix A 



All the splitting functions in the Altarelli-Parisi (AP) equations can be ex- 
panded as a perturbation series in ag into LO and NLO terms as follows 

P^p(o) + ^p(i). (A.i) 
27r 

The non-singlet combinations of the qr{qr) to q's(^s) splitting functions, where 
the subscripts r, s denote the flavors of the (anti)quarks and satisfies r, s = 
1, ■ ■ ■ , n/, can be further decomposed into a flavor diagonal part proportional 
to 5rs a-nd a flavor independent part. In LO there is only one non-singlet 
splitting function Pqq but in NLO it is convenient to form two combinations 
from Pqq and Pqq as follows 

P+ — Pqq + Pqq 
P_ = Pqq - Pqq. (A.2) 

These splitting functions are used to evolve two independent types of non- 
singlet densities, which will be called plus and minus respectively. They are 
given by 

f~ = fk^k{nf,x,fj,^). (A.3) 

The easiest way to explain the indices in these equations is by explicitly 
giving the combinations we use. For j = 1,2 

f-^u-u,f2=d-d, (A.4) 

which are used for all flavour density sets. Then for three-flavor densities 
i — 1,2,3 and we define 

f+^u + u- E(3)/3 , f+^d + d- E(3)/3 , 

/H-^s + ^-E(3)/3, (A.5) 

where S(3) = /^(3) = u + u + d + d + s + s. These densities should be used 
for II < nic- For four-flavor densities i — 1, 2, 3, 4 and we deflne 

f+ = u + u- S(4)/4 , f+ = d + d- S(4)/4 , 

f+ = s + s- S(4)/4 , /+ = c + c - S(4)/4 , (A.6) 
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where E(4) = /^(4) = c + c + S(3). These should be used for rUc < A* < rrib. 
For five-flavor densitcs i = 1, 2, 3, 4, 5 and we define 



/+ =-ii + -u-E(5)/5, /2+ = (i + d- E(5)/5, 
/3+ = s + 5-S(5)/5, /4+ = c + c-E(5)/5, 
/5+ = 6 + 6-E(5)/5, 



(A.7) 



where E(5) = /|(5) = 6 + ft + E(4). These should be used for // > m^. 

If we define t = ln(/x^/(l(GeV/c^)^) then the AP equations that we need 
to solve are 



dt 



(A.8) 



dt 



(A.9) 



dfgix,t) 
dt 



as{t) dz 



I 

Jx 



27r Jx z 



pA^i^,t) + p,Af^{z,t) 



, (A.IO) 



df^(x,t) _ as{t) f^dz 



dt 



/ 

Jx 



27r Jx z 



PU-)f,'i^,t) + P,,{-)f^{z,t) 

z z . 



, (A.ll) 



where for /j, < rUc we set i — 1,2,3, j — 1,2, — E(3) and the gluon is 
a three-fiavor gluon. When rric < < mi,, we use i = 1,2,3,4, j — 1,2, 

= E(4) and the gluon is a four-fiavor gluon. Finally when > m^, we set 
i = 1,2. 3, 4,5, j = 1, 2, fji = E(5) and the gluon is a fivc-fiavor gluon. 

The densities should satisfy the momentum conservation sum rule 



/ X u(x, ii"^) + d(x, 11^) + u(x, 11^) + d(x, ii'^) + 2s(x, 11^) 
Jo L 



+2c{x, n^)9{ii^ - m^) + 2b{x, n^)9{fi^ - ml) + g{x, fi^) 



dx = 1 . 

(A.12) 



As the quark constituents carry all the charge, isospin, strange, charm and 
bottom quantum numbers of the nucleon they should also satisfy the other 
standard sum rules for the conservation of these quantities. 
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There are several methods to solve these equations. Among them the 

see full list of 



23 , 24 



most popular are to use Mellin moments (used by 
references in [|IB| and [|T^) and to use the direct x-space solution (as in 
21 



17 



see also [|T8| ). Also the authors in [^] describe a method involving 
Laguerre polynomials, that dates back to early paper of 0. 

Our choice of direct x-space method is justified by the necessity to step 
densities across matching points using LO, NLO and NNLO boundary con- 
ditions. The procedure of doing this in the Mellin moment method would 
involve converting densities to and from Mellin moments several times. Us- 
ing the direct x-space method is much more intuitive and straightforward. 
The main features of this method are linear interpolation over a grid in x 
and second-order interpolation over a grid in t. Let us describe the method 



in more detail to point out where we differ from the work in ||T7|. 

First we consider the x-variable in the evolution. Consider the right- 
hand-side of the evolution equation ( [A.8| ) for non-singlet density 



/(xo) 



dz Xq 



z z 



p 



Xq 



q{z) 



where xq < z <1 and 
and 



q[x) = xf{x) , 



(A.13) 



(A. 14) 



Xq<Xi< ... <Xn< Xn+l = 1 , (A.15) 

with q{xn+i) = qiX) = 0. Between grid points Xi and Xj+i, x is chosen so 
that 

q(^x) = {1 - y)q{x,) + yq{xi+i) , (A.16) 

with y = (x — Xi)/ (xj+i — Xj). Using this relation we convert the integral into 
a sum 



n+l 



H^o) = w{xi,Xo)q{xi) , 

i=0 



(A.17) 



where the weights are 



w{xo,xo) = S'i(si,So) 

w{Xi,Xo) = Si{Si+i,Si) - S2{Si,Si^i] 



(A.18) 
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where Si = xq/xi and 

Si{u,v) 
S2{u,v) 



V 



V — u 
u 

V — U Ju 



V rj'/ 

{z-u)P{z) — 

{z-v)P{z)—. 

z 



(A.19) 



We have calculated these integrals analytically and the results are in the 
computer program. This leads to the final formula describing the grid for 
the X variable. Note that the weights and w^^^ include LO and NLO 
splitting functions respectively. Thus, for the singlet case, we have 



(i(xoS(xo)) as 



dt 



2ti 

+ 



E 



wf^i.Xi,XQ) + ^wW(Xi,Xo) 



Xig{xi) . (A.20) 



Now consider the variation in the variable t. For each Xj we pick a grid in 
t labelled by distinct points tj. Then the example the non-singlet equation 
becomes 



k=l 



W^^\xk,Xi) + 



2ti 



q{xk,t,), (A.21) 



where q{xi,tj) denotes the derivative with respect to t evaluated at t = tj. 
In compact notation this equation can be rewritten as 



q'j = wqj + S , 



(A.22) 



with S being the sum of the terms on the right hand side of ( [A. 211) excluding 
the j-th term. 

For t between the grid points and tj we interpolate the parton density 
using quadratic interpolation as follows 

q{x^,t) = at^ + bt + c. (A.23) 

Thus we relate the value of q at the point tj to that of q at the point tj^i by 

q{xi,tj) = q{xi,tj-i) + ^[q {xi,tj) + q {xi,tj^i)]Atj , (A.24) 
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where Atj = tj — tj-i. This equation can also be written more compactly as 



1 



= U-i + 2^<lj-i + (Ij)^t, . (A.25) 



The resulting system of two linear equations ( |A.25|) and (|A.22|) for qj and 
has the solution 

and yields from ( |A.22|) . Applying the same procedure to the gluon and 
singlet equations Eqs.(A.10)-(A.ll) involves four equations because we have 
to compute both the densities and their derivatives. 

The evolution proceeds from the initial /ig = /Ulo (^^ lA — /^nlo) 
the first matching point at the scale //^ = m^. Next the charm density is 
introduced in the NNLO (a^-order terms) and all the four-flavor densities 
are evolved from the boundary conditions in Ens. (|2.8| )- (|2.1CI| ). This evolution 
continues up to the transition point /i^ = m^, where the same procedure 
is applied to generate the bottom quark density. At that matching point 
all five-flavor densities are evolved starting from the boundary conditions in 
Eqs.( p.8|) -( |2.11|) up to all higher scales. 
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Figure Captions 

Fig. 1. (a) The charm quark density xcnnlo(4, x, jji^) the range 10~^ < x < 
1 for = 1.96, 2, 3, 4, 5, 10 and 20 in units of (GeV/c^)2, (b) similar 
plot as in (a) but now for 0.01 < a; < 1, (c) ratios yU^) = 

2^cevolved(4, a;, /x^)/a;cFOPT(4, a;, /x^) for the same scales, (d) and (e) 
the NLO results from MRST98 set 1 and CTEQ5HQ respectively. 

Fig. 2. (a) The gluon density a;(7NNLo(4, /i^) in the range 10^"'^ < x < 
1 for /i^ = 2, 3, 4, 5, 10 and 20 in units of (GeV/c^)2, (b) ratios 
^NNLO^^^ /i^) = a;(7EVOLVED(4, X, /x^)/a;(7FOPT(4, x, iJ?) for the same scales, 
(c) the three-flavor NLO gluon density in the same range, (d) and (e) 
the NLO results from MRST98 set 1 and CTEQ5HQ respectively. 

Fig. 3. (a) The singlet density a;ENNLo(4, a:;, /^^) in the range 10~^ < x < 
1 for /i^ = 2, 3, 4, 5, 10 and 20 in units of (GeV/c^)^ (b) ratios 
j^NNLOj^^^ = a;EEvoLVED(4, X, /a;EFOPT(4, x, //^) for the same 
scales, (c) the three-flavor NLO density. 

Fig. 4. (a) The nonsinglet quark density X(Tnnlo(4, x, ij?) in the range 10~^ < 
a; < 1 for = 2, 3, 4, 5, 10 and 20 in units of (GeV/c^)2, (b) ra- 
tios i?^^^^(a;, /i^) = X(Tevolved(4, X, /i^)/a;(jFOPT(4, a;, /x^) for the same 
scales, (c) the three-flavor NLO density. 

Fig. 5. (a) The strange quark density a;sNNLo(4, a;, /i^) in the range 10~^ < 
a; < 1 for /x^ = 2, 3, 4, 5, 10 and 20 in units of (GeV/c^)^, (b) ra- 
tios /x^) = a;cEvoLVED(4, a;, /x^)/a;cFOPT(4, a;, /i^) for the same 
scales, (c) the three-flavor NLO density. 

Fig. 6. (a) The bottom quark density a;6NNLo(5, a;, /i^) in the range 10^^ < 
a; < 1 for = 20.25, 25, 30, 40 and 100 in units of {GeV/cf, (b) sim- 
flar plot as in (a) but now for 0.01 < a; < 1, (c) ratios B^^^^{x,n'^) = 
x6evolved(4, a;, /i^)/x6FOPT(4, x, ji^) for the same scales, (d) and (e) the 
NLO results from MRST98 set 1 and CTEQ5HQ respectively 

Fig. 7. (a) The charm quark density a;cNNLo(5, a;, //^) in the range 10~^ < 
a; < 1 for = 20.25, 25, 30, 40 and 100 in units of (GcV/c^)2, (b) and 
(c) the NLO results from MRST98 set 1 and CTEQ5HQ respectively 
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Fig. 8. (a) The gluon density xg-^^i^o{5, x, fi"^) in the range 10 < a; < 1 
for /i^ = 20.25, 25, 30, 40 and 100 in units of (GeV/c^)^, (b) the three- 
flavor NLO density. 

Fig. 9. (a) The singlet quark density a;SNNLo(5, a;, /x^) in the range 10~^ < 
x < 1 for = 20.25, 25, 30, 40 and 100 in units of (GeV/c^)^ (b) the 
three-flavor NLO density. 

Fig. 10. (a) The nonsinglet density a;(jNNLo(5, a;, /x^) in the range 10~^ < 
X < 1 for ^ 20.25, 25, 30, 40 and 100 in units of (GeV/c^)^ (b) the 
three-flavor NLO density. 

Fig. 11. (a) The strange quark density a;sNNLo(5, x, fi^) in the range 10^^ < 
X < 1 for ^ 20.25, 30, 40 and 100 in units of (GeV/c^)^, (b) the 
three-flavor NLO density. 
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